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TITLE OF THE INVENTION 
System and Method For Diagnosing Pathologic Heart Conditions 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims the benefit of prior filed co-pending U.S. Provisional Patent 
Application No. 60/457,726, filed on 3/26/03. 

BACKGROUND OF THE INVENTION 
[0002] The present invention relates to a system and method for diagnosing pathologic heart 
conditions based upon heart sound data. 

[0003] Studies have shown that primary care physicians frequently refer patients to cardiac 
specialists on the basis of suspicious heart sounds detected by traditional stethoscope 
auscultation, though a large percentage of these referrals are dismissed by cardiologists as having 
no pathologic condition. The costs, delays, worry, and administrative burden resulting from 
these needless referrals could be reduced if the cues that the specialist uses could be incorporated 
into an algorithm to automatically screen for pathologic heart sounds and murmurs. 
[0004] Although attempts have been made to automate screening by auscultation, no device is 
currently available to fulfill this function. Multiple indicators of pathology are nonetheless 
available from heart sounds and can be elicited using certain signal processing techniques such 
as time-frequency analysis. Some pathology signals can reliably be detected and classified as 
pathologic using a portable electrocardiogram and heart sound measurement unit combined with 
time-frequency based algorithms. Time-frequency decomposition analysis holds promise for 
extending these results to detection and evaluation of other audible pathologic indicators. 
[0005] In addition, an automatic screening algorithm would be useful for detecting pathologic 
heart conditions in settings where a trained professional is not always available, such as pre 
sports participation physicals, and examinations performed in remote or underserved areas. 
Furthermore, automated analysis of digitized clinical information such as heart sounds could 
have major implications for health care delivery systems using telemedicine. 
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SUMMARY OF THE INVENTION 
[0006] The present invention comprises a time-frequency diagnostic device and method. The 
present invention integrates a cardiologist's auscultation expertise, a large and growing set of 
comprehensive heart sound files, and digital signal processing algorithms. 
[0007] A method of diagnosing pathologic heart conditions in which a time series of heart 
sounds and ECG peaks is filtered and parsed into a sequence of individual heart cycles. A 
systolic interval as well as systolic sub-intervals are identified for each heart cycle. 
[0008] The systolic intervals are then digitally filtered to optimize for click detection. A digital 
filter transform output is then fed to multiple processes that occur for each heartcycle including 
determining the systole time limits, inputting the time series of the transformed heartsound 
output to a Neyman-Pearson "constant false alarm rate" (CFAR) detector to identify 
anomalously high transform coefficients, and creating a vector of detections vs. time. The series 
of anomalously high detections (one series for each heart cycle) are then assembled into a matrix. 
Each heart cycle in the matrix is convolved with an averaging vector whose averaging length is 
consistent with an observed spread of click occurrence times in a database of patients. . A click 
score is then determined from these averaged detections as the maximum element of the vector 
formed by the median transform coefficient amplitude across heart cycles squared at each time 
sample multiplied element-wise by the vector formed by the sum across heart cycles of the 
number of detections at each time sample. 

[0009] The click score is compared to a threshold value set by a desired probability of detection 
vs. a probability of false alarm tradeoff. If the click score is less than the threshold then a "no 
click" indicator is displayed. If the click score is greater than the threshold then a "click present" 
indicator is displayed. The click score of the patient in question is compared to the distribution of 
scores of known pathologic cases and known non-pathologic cases from the database. The 
ranking of the current score in the distribution of pathologic and non-pathologic scores is 
presented to the operator for additional interpretation context . 

[0010] The system for diagnosing pathologic heart conditions is comprised of a portable 
computing device that manages data collection and stores data collected from new patients, and 
analyzes data. Also included is a patient data collection unit, communicable with the portable 
computing device, for acquiring and digitizing electro-cardiogram (ECG) and heart sound data 
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from a patient. The patient data collection unit is comprised of a pair of transducer contact 
microphones (primary and reference) for obtaining acoustic data. Also included is a pair of 
acoustic pre-amplifiers connected with the transducers. The pre-amplifiers have a passband of 20 
Hz to 2 kHz and are used to condition acoustic data received from the contact microphones. 
Variable amplifiers connected with the acoustic pre-amplifiers variably amplify the conditioned 
acoustic data. Moreover, several electro-cardiogram (ECG) electrodes connected to an ECG 
amplifier record ECG data. The acoustic and the ECG data are passed to an analog to digital 
converter connected with the variable amplifiers and the ECG amplifier. The data is digitized 
and sent to the computing device for processing by a screening algorithm implemented by the 
steps described in the method above. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0011] FIGURE 1 illustrates a plot of time-frequency analysis at various scales a. 
[0012] FIGURE 2 illustrates a block diagram of system hardware. 

[0013] FIGURE 3 illustrates a logic flow diagram of the processes used to diagnose pathologic 
heart conditions such as heart clicks. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0014] Before describing the present invention it is helpful to have a basis for the detection of 
pathological heart conditions. Analysis of a variety of heart sounds and corresponding diagnoses 
from the Johns Hopkins School of Medicine (JHU/SOM) has shown that time-frequency analysis 
is a versatile technique for detecting and classifying pathologic heart conditions. Of the 
available time-frequency techniques, wavelets are a useful method for representing heart sound 
frequency dynamics without creating cross term artifacts. Wavelet transforms can be computed 
for a continuous or discrete set of transform variables, depending on the priority for compactness 
(e.g., data compression applications) or ease of visual interpretation (e.g., pattern recognition 
applications), respectively. While compressibility is of interest for future applications, the 
relative ease of interpretation provided by continuous wavelet transforms (C WTs) suggested 
their use with the present invention. While wavelet transform analysis are illustrated herein as a 
preferred method for analyzing heart sound data, frequency based Fourier transform analysis 
may also be implemented by one of ordinary skill in the art in order to analyze heart sound data. 
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where f and v|/ are both square-integrable, a is a time scaling variable, and b is a time translation 
variable. This can also be written as a convolution: 



[0016] While it is possible to construct a v|/ to yield an optimal peak (compact, high amplitude) 
in W for a given f, it would not be guaranteed to be optimal for a different time series g. Since 
the present invention is to be applied to a variety of heart sound signals, a custom wavelet was 
not implemented. Rather, the alternative was to draw vj; from a pool of wavelets designed to 
have various advantageous properties. Wavelet transformations known as second order 
"coiflets" were implemented with the present invention. It is important to note, however, that 
other wavelet transformations, including custom wavelets, may be implemented without 
departing from the spirit or scope of the present invention. 

[0017] Coiflets of order 2 at various scales, a, are plotted in FIGURE 1. The frequency 
bandpass limits at each scale are given. Due to the equivalence of convolution in the time 
domain and multiplication in the frequency domain equation 1 shows that the Fourier transform 
of \]f will act as a bandpass filter of the signal f. The bandpass limits in FIGURE 1 are the 6 dB 
passband frequency limits of the Fourier transform of (p . Wavelets are constructed so as to 
maintain a constant ratio of center frequency to 3 dB bandwidth (Q), and have a finite duration. 
Their time-frequency resolution is inherent in their design and scale parameters. This is in 
contrast to Fourier decomposition, which uses the infinite time extent sine and cosine functions. 
Time resolution is not inherent in the Fourier transform, but is introduced by the user via 
windowing the data. Multiple Fourier transforms using distinct window intervals would be 
required to produce the constant Q decomposition offered by wavelets. 

[0018] A system block diagram of the present invention is illustrated in FIGURE 2. The system 
is comprised of two principal elements. One is a patient data collection unit 201. The other is a 
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computer processing device 220 including or having access to data storage devices. The patient 
data collection unit performs several functions including obtaining heart sound data via a set of 
contact microphones in the form of transducers. Two channels of acoustic data are obtained 
from a patient using a primary transducer 202 and a reference transducer 204. In addition, a set 
of ECG electrodes 206 are used to obtain electro-cardiogram data from the patient. 
[0019] The two contact microphones are each conditioned by a pre-amplifier 208 having a 
passband of 20 Hz to 2 kHz and variable gain amplification stage 210. A set of headphones 212 
connected to a headphone amplifier 214 can be used to listen to the acoustic data gathered from 
the primary transducer 202 and the reference transducer 204. The ECG electrodes 206 feed into 
an ECG amplifier 216. Outputs from the variable gain amplifiers 210 and the ECG amplifier 
216 are fed to a analog-to-digital converter where the acoustic and ECG signals are digitized and 
recorded. A 12-bit National Instruments PCMCIA analog-to-digital converter is used, for 
instance, to collect and digitize data at a rate (e.g. 8.13 kHz) consistent with the highest data 
frequencies of interest. 

[0020] Once the signals have been digitized and recorded, the patient data collection unit 201 
forwards the data to a computer processing device 220. The computer processing device 220 is 
typically a laptop computer (due to its compact transportable nature) having adequate data 
storage capacity. However, the patient data collection unit 201 may be connected to other 
computer processing devices without departing from the spirit or scope of the present invention. 
A computer software program accesses heart sound data that has either been collected and 
forwarded by the patient data collection unit 201, or is resident on the laptop computer 220, or 
can be obtained from another source of heart sound data. 

[0021] The computer program applies a screening algorithm to the heart sound data in order to 
determine whether the heart sound data is to be classified as normal or pathologic. Each heart 
sound data file corresponds to a different patient. When the algorithm has operated on the heart 
sound data that has been input, the computer program will display the results on a display screen 
to the doctor, nurse, or technician operating the computer. Results indicating a pathologic 
condition will likely cause the patient to be referred to a cardiologist for further examination. 
Otherwise, a cardiologist referral can be deemed unnecessary. 

[0022] FIGURE 3 illustrates the flow of logic and processing that occurs in the various elements 
described in FIGURE 2. Digitized heart sound recordings were collected on patients in the 
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Pediatric Cardiology Echocardiography Laboratory in the Johns Hopkins Outpatient Center. 
The recordings were stored in a Heart Sound database 302. Recorded heart sounds could then be 
extracted 304 from the heart sound database 302 and placed into an ECG/heart sound data file set 
306. Alternatively, ECG and heart sound data could be obtained directly from a patient 308 and 
placed into the ECG/heart sound data file set 306. The ECG/heart sound data file set 306 serves 
as the data to be fed to a screening algorithm. The purpose of the screening algorithm is to 
analyze the ECG data and heart sound data in order to detect pathologic anomalies that may be 
present. Thus the system is to be used as a diagnostic aid. In order for the screening algorithm 
to be applied, the data set must first be manipulated. Initially, the ECG data and acoustic heart 
sound data are separated. The ECG data is used to parse the time series into a sequence of 
individual heart cycles via a process that identifies the ECG peaks 310. 

[0023] A systolic interval is then identified 312 for each heart cycle. The first and second heart 
sounds are identified both by reference to the heart cycle boundaries and acoustically using a 
passband of 25-140 Hz 314. In the acoustic method the times of the acoustic maxima define 
systole and diastole. Systole can then be divided into various sub-intervals including, but not 
limited to, the first half of the systolic interval, the total systolic interval, the middle half of the 
systolic interval, and the last half of the systolic interval. A short interval to isolate the first and 
second heart sounds is factored in. Then, the acoustic heart sound data is passed through a 
digital bandpass filter 316, in this case a second order coiflet CWT transform. The purpose of 
this filter is to optimize the acoustic data series for click detection. The digital wavelet filter 
output of box 316 is then fed to multiple processes that occur for each heartcycle 318. 
[0024] For each heartcycle 318, the systole time limits are determined, a time series of the 
transform at specific wavelet scales (16 and 24 for instance) are input to a Neyman-Pearson 
"constant false alarm rate" (CFAR) detector to identify anomalously high wavelet coefficients, 
and a vector of detections vs. time is created. The noise model for the Neyman-Pearson CFAR 
detector is statistically based on the study of heart sound wavelet coefficients in a variety of 
patients. It was determined to be best fit by a gamma distribution. 
[0025] The series of anomalously high detections (one series for each heart cycle) are then 
assembled into a matrix 320. The matrix is convolved with an averaging vector, whose length is 
derived from the expected time spread of click occurrence yielding detection statistics across 
heart cycles and time intervals consistent with an observed spread of click occurrence times 322. 
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A click score is then calculated 324. The click score is determined as the maximum element of 
the vector formed by the median wavelet coefficient amplitude across heart cycles squared at 
each time sample multiplied element-wise by the vector formed by the sum across heart cycles of 
the number of detections at each time sample. 

[0026] The click score is fed to a decision box 326 in which it is compared to a threshold value. 
The threshold value is set by a desired probability of detection vs. a probability of false alarm 
tradeoff. If the click score is less than the threshold then a "no click" indicator 328 is displayed 
334. If the click score is greater than the threshold then a "click present" indicator 330 is 
displayed 334. 

[0027] In addition, the click score is compared to other click scores in a database and the newly 
calculated click score is presented 332 relative to other click scores that are both normal and 
pathologic. The newly calculated click score and all the corresponding data that led to its 
calculation is saved in a patient database 336. The patient database continues to grow and 
provide greater reliability and accuracy as it becomes populated with more and more samples. 
[0028] As the data is being processed certain conditions are monitored 338, 340 by the present 
invention to indicate whether there is an error in processing. The error conditions presently 
monitored include: a) that the ecg data allows parsing of the heartcycles; b) reasonable values are 
obtained for SI and S2 peak location; and c) that the length of systole is enough to contain the 
split window used for the CFAR processor. 

[0029] There are several advantages realized by the present invention. An automatic screening 
algorithm would be useful for detecting pathologic heart murmurs in settings where a trained 
professional is not always available, such as pre sports participation physicals, and examinations 
performed in remote or underserved areas. Furthermore, automated analysis of digitized clinical 
information such as heart sounds could have major implications for health care delivery systems 
using telemedicine. 

[0030] Automated analysis of heart sound data could also be used by cardiologists to 
quantitatively follow and document changes in severity of certain conditions such as aortic 
stenosis and mitral regurgitation, in which changes in murmur characteristics are known to 
correlate with changes in disease severity. In addition, since data for analysis could theoretically 
be collected using any electronic stethoscope, home care nurses or nurses aides could make 
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inexpensive in-home bedside recordings that could be analyzed later to monitor changes in 
certain conditions as reflected in heart sounds. 

[0031] It is to be understood that the present invention illustrated herein is readily implementable 
by those of ordinary skill in the art as a computer program product having a medium with 
computer program(s) embodied thereon. The computer program product is capable of being 
loaded and executed on the appropriate computer processing device(s) in order to carry out the 
method or process steps described. Appropriate computer program code in combination with 
hardware implements many of the elements of the present invention. This computer code is 
typically stored on removable storage media. This removable storage media includes, but is not 
limited to, a diskette, standard CD, pocket CD, zip disk, or mini zip disk. Additionally, the 
computer program code can be transferred to the appropriate hardware over some type of data 
network. 

[0032] The present invention has been described, in part, with reference to flowcharts or logic 
flow diagrams. It will be understood that each block of the flowchart diagrams or logic flow 
diagrams, and combinations of blocks in the flowchart diagrams or logic flow diagrams, can be 
implemented by computer program instructions. 

[0033] These computer program instructions may be loaded onto a general purpose computer, 
special purpose computer, or other programmable data processing apparatus to produce a 
machine, such that the instructions which execute on the computer or other programmable data 
processing apparatus create means for implementing the functions specified in the flowchart 
block or blocks or logic flow diagrams. 

[0034] These computer program instructions may also be stored in a computer-readable memory 
that can direct a computer or other programmable data processing apparatus to function in a 
particular manner, such that the instructions stored in the computer-readable memory produce an 
article of manufacture including instruction means which implement the function specified in the 
flowchart blocks or logic flow diagrams. The computer program instructions may also be 
loaded onto a computer or other programmable data processing apparatus to cause a series of 
operational steps to be performed on the computer or other programmable apparatus to produce a 
computer implemented process such that the instructions which execute on the computer or other 
programmable apparatus provide steps for implementing the functions specified in the flowchart 
blocks or logic flow diagrams. 
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[0035] Accordingly, block(s) of flowchart diagrams and/or logic flow diagrams support 
combinations of means for performing the specified functions, combinations of steps for 
performing the specified functions and program instruction means for performing the specified 
functions. It will also be understood that each block of flowchart diagrams and/or logic flow 
diagrams, and combinations of blocks in flowchart diagrams and/or logic flow diagrams can be 
implemented by special purpose hardware-based computer systems that perform the specified 
functions or steps, or combinations of special purpose hardware and computer instructions. 
[0036] In the following claims, any means-plus-function clauses are intended to cover the 
structures described herein as performing the recited function and not only structural equivalents 
but also equivalent structures. Therefore, it is to be understood that the foregoing is illustrative 
of the present invention and is not to be construed as limited to the specific embodiments 
disclosed, and that modifications to the disclosed embodiments, as well as other embodiments, 
are intended to be included within the scope of the appended claims. The invention is defined 
by the following claims, with equivalents of the claims to be included therein. 
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